Is there evidence for sterile neutrinos in IceCube data? 
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Data from the LSND and MiniBooNE experiments, and revised expectations of the antineutrino 
flux from nuclear reactors suggest the existence of eV-mass sterile neutrinos. 3 + 2 and 1 + 3 + 1 
scenarios accommodate all relevant short-baseline neutrino data except for the low-energy Mini- 
BooNE anomaly. We analyze the angular distribution of upward going atmospheric neutrino events 
in the IceCube-40 dataset for evidence of sterile neutrinos within these scenarios. Depending on how 
systematic uncertainties are handled, we find strong evidence for, or weak evidence against sterile 
neutrinos. We show that future IceCube data will definitively settle the issue. 



While oscillations of the three Standard Model neutri- 
nos have been convincingly demonstrated, and the frame- 
work that explains the data of many disparate experi- 
ments is elegant in its simplicity, a nagging doubt per- 
sists. Are we seeing the effects of light sterile neutri- 
nos in short-baseline neutrino (SBL) experiments? The 
question is raised by 3.3cr CL evidence from the Liquid 
Scintillator Neutrino Detector (LSND) [I] experiment for 
f'/i ^ i^e oscillations with L/E ^ 1 ni/MeV. This result is 
partially corroborated by the Mini-Booster Neutrino Ex- 
periment (MiniBooNE) [2]. The mass-squared difference 
6m^ ^ 1 eV^ indicated is very different from those that 
explain solar and atmospheric neutrino data. Since three 
neutrinos have only two independent mass-squared dif- 
ference scales, a third 6m^ scale implies the existence of 
additional light neutrinos, which must be sterile because 
the invisible width of the Z boson requires that there be 
no more than three active neutrinos. 

A new theoretical calculation of the J>e flux from nu- 
clear reactors finds a value that is 3% higher than previ- 
ous estimates [3]. If borne out, the incident flux at SBL 
reactor neutrino experiments may be interpreted as hav- 
ing a deficit due to oscillations into eV sterile neutrinos. 

A recent analysis of SBL beam and reactor experi- 
ments, using the new reactor flux prediction, concludes 
that a 3 + 1 scenario with a single additional eV sterile 
neutrino does not describe the data satisfactorily, while 
a 3 + 2 scenario with two eV sterile neutrinos, and a 
1 + 3 + 1 scenario with one sterile neutrino lighter than 
the 3 active neutrinos and the other heavier, provide a 
good fit; see Table |l] for the best- fit parameters [3]. The 
phase S = arg[U^iU*/JJ'^n,Ue5), where U is the neutrino 
mixing matrix, leads to CP violation. Note that accom- 
modating such eV sterile neutrinos in cosmology requires 
quite drastic modifications of ACDM (5j . 

The IceCube experiment (IC) [S] observes lO'^ atmo- 
spheric muon neutrino events per year thereby offering 
a unique probe of neutrino oscillations at energies above 
100 GeV. For upward going neutrinos, any additional 
species that significantly mixes with impacts the ob- 
served muon rate through disappearance. In fact, for 



(?mli 5mli \Ue4\ \UtJ,A\ \Ue5\ \U^ib\ 

3+2 0.47 0.87 0.128 0.165 0.138 0.148 1.64 

l+3+l(a) -0.47 0.87 0.129 0.154 0.142 0.163 0.35 

l+3+l(b)| 0.47 -0.87 I 0.129 0.154 0.142 0.163 1 1.65 

TABLE I. Global best-fit parameters to data from short- 
baseline experiments 4-. The two 1 + 3+1 cases correspond to 
either m4 or ms being the lightest state. The active neutrinos 
have a normal hierarchy, and ^13 = 0. 

5m? ~ 1 eV^, resonant oscillations are possible for TeV 
atmospheric neutrinos [71. In its 40-string configuration, 
IC has observed 12877 atmospheric muon neutrino events 
in the energy range 332 GeV-84 TeV 8J. 

In this Letter, we analyze the zenith-angle distribution 
of the IC data to see if it provides supporting evidence 
for sterile neutrino explanations of the SBL data and the 
reactor antineutrino anomaly as encapsulated by the sce- 
narios of Table |Tj and assess the future sensitivity of the 
IC detector. The energy window of the IC dataset is op- 
timal as can be seen from Fig. [T] We expect IC data 
to not be very sensitive to the 3 + 2 case since the reso- 
nance occurs in the channel which is subdominant in 
the atmospheric flux. In the 1 + 3 + 1 cases, the domi- 
nant i/^ flux is also suppressed so that these cases can be 
discriminated more readily from the 'Hv case. 

The muon events at IC can be classified as (i) 'con- 
tained' events, namely events with muon tracks that 
start within the instrumented volume, and (ii) 'up-going' 
events in which the muon is produced outside the detec- 
tor. See Eqs. 8 and 10 of Ref. [lOj for details and note 
that the incident i^^ flux includes a contribution from 
i^e — >■ oscillations. 

For the points in Tablejlj the muon neutrino flux is sup- 
pressed by ~ 10% above 100 GeV, with a corresponding 
reduction in the muon event rate. Thus, the 25% uncer- 
tainty in the atmospheric neutrino flux |llj is a serious 
impediment to a rate analysis. Instead, a distortion in 
the zenith-angle distribution of the detected muons would 
provide the strongest evidence for sterile neutrinos. As 
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FIG. 1. The survival probability for j/^ (upper panel) 
and (lower panel) for vertically going upward neutrinos 
(cos^z = — 1) in the 3 scenarios of Table |l] We use the den- 
sity profile of the Preliminary Reference Earth Model [9] . 



most sub-TeV neutrinos do not distort the angular distri- 
bution, useful information may be extracted from high- 
energy events. It so happens that the cut at 332 GeV 
suppresses 90% of sub-TeV muon events while leaving a 
large sample of 10"* events. 

To make it to the detector, up-going muons can travel 
at most the stopping distance which increases with en- 
ergy. As a result, the effective volume for neutrinos is 
larger for high energy neutrinos that more readily pro- 
duce energetic muons. In turn, the up-going sample has 
a larger fraction of high-energy events than the contained 
sample, making it particularly valuable in an analysis of 
the zenith-angle distribution. The 40-string IC data con- 
sist of 70% up-going and 30% contained events [12]. 

To account for experimental efficiencies, we determine 
the theoretical zenith-angle spectrum as in Ref. [T5] : 



5f = a[l + fe((cos0,),+O.5)]^7Y. 



th 



(1) 



where Nf^ is the number of muon events in angular bin 
i, and N^^'-^ /Nf'^ is a bin- wise factor that scales our the- 
oretical iv prediction to IceCube's Monte Carlo [5]. a is 
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FIG. 2. Upper panel: Zenith-angle distributions for muon 
events at IC. The IC data (light, orange) and the expectation 
for 2>v oscillations (dark, blue) are shown with 2o statistical 
uncertainty bands. The expectations for the 3 + 2 (dark, 
red), 1 -I- 3 -I- 1(a) (light, green sohd), and 1-f 3 -H(b) (light, 
green dotted) scenarios are shown. We set a = 1 , & = to 
display the rate suppression and the spectral shape in each 
case. Lower panel: Relative distortion of the zenith-angle 
spectra in comparison to the expectation for oscillations. 
Here, the best-fit value of a is used in each case, but 6 = 0. 
The data reveal an obvious distortion with respect to the ?>u 
expectation. 



an overall normalization and b allows a systematic tilt of 
the spectrum. Both parameters are allowed to float in a 
2 analysis. Note that Sf" = TV*^"^ in the iv case with 



1,6 = 0. 



The upper panel of Fig. [2] illustrates the event rate 
suppression and zenith-angle distributions in the scenar- 
ios under consideration. The large difference of the sterile 
neutrino scenarios from the IC data and from the 2>v ex- 
pectation in the near-vertical (—1.0 < cos^z < —0.9) bin 
promises much discriminating power. The near-horizon 
(—0.1 < cos < 0) bin has potentially large systematic 
errors from the misidentification of coincident downward 
events as horizontal events, and we do not include it in 
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FIG. 3. x'^ Eis a function of the tilt parameter b from a fit to 
the 9 zenith-angle bins of Fig. [2] with cos^^ < —0.1. 



our analysis. For our analyses we define 



(1 



th 



(2) 



where <7a = 0.25 is the percent uncertainty in the atmo- 
spheric flux normalization [TT] and S^^^ denotes either 
real or simulated data. In what follows, we always fit a, 
and either set & = or marginalize over b without penalty. 
The lower panel of Fig. [2] shows the spectral distortion of 
the IC data and the sterile neutrino scenarios relative to 
Sf oscillations. In the latter case, a = 1 , 6 = 0, while the 
best-fit value of a is used for the sterile cases with 6 = 0. 

The results of fitting the 9 non-horizontal (cos^z < 
—0.1) bins are shown in Fig. [S] and Table [lT| The sce- 
nario gives a better fit provided the data are plagued by a 
large systematic tilt. If |&| is restricted to be smaller than 
0.11, the Siy fit yields a of 15.8, which is comparable 
to the values (~ 15-16) of the sterile neutrino cases. 

A natural question is if there are correlated signals 
in cascade events at IC. Interestingly, above 332 GeV, 
we find a comparable total number of Ve plus Vr 
events/km^ /year in the Siy 890) and sterile neutrino 
(~ 850) cases. The corresponding number of induced 
cascade events can be obtained by folding with the IC 
efficiency which is not available to us. 

That future IC data with systematics under control 
have the ability to reveal a deviation due to sterile neu- 
trinos is evident from Fig. |4j The monotonically rising 
event ratio as a function of cos 9z in the up-going event 
sample is a striking signature of sterile neutrino oscil- 
lations. On the other hand, strong exclusions can also 
be obtained. Assuming a sample of 6.5 x lO'' up-going 
events with no deviation from the 3^ result, and fitting 
to all 10 angular bins with only a varied yields = Hlj 
386 and 331 for the 3-1-2, 1 -I- 3 -I- 1(a) and (b) cases, 
respectively, compared to a x'^ ^ 10 for oscillations. 
An interesting aspect of the contained event sample is 
that in all cases the zenith-angle distribution is almost 
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35.2 


13.0 


0.17 


3 + 2 


18.5 


15.2 


0.066 




18.1 


15.6 


-0.058 


1 + 3+ 1(b) 


16.4 


16.0 


-0.024 



TABLE II. Best-fit and b from the curves in Fig.|3] With- 
out a systematic tilt in the angular distribution, sterile neu- 
trino scenarios are clearly favored over the case. If a large 
systematic tilt is permissible, the sterile cases are mildly dis- 
favored. For |6| < 0.11, all scenarios give a comparable fit. 
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FIG. 4. Similar to the lower panel of Fig. [2j but with 
1.3 X 10"" events above 332 GeV equipartitioned into purely 
up-going events (left panel) and purely contained events (right 
panel). The shaded band represents the 2a statistical uncer- 
tainty on the expectation. Sterile neutrino scenarios will be 
easily distinguishable from oscillations from the up-going 
event sample. 



flat for —0.8 < cos 6*2 < —0.2 so that the features in the 
near-vertical and near-horizon bins are insensitive to a 
systematic tilt. With 6.5 x 10^ contained events, we find 
X^ = 28, 165 and 87 for the 3 + 2 and 1 + 3 + 1(a) and (b) 
cases, respectively. With such sensitivity IC will easily 
confirm or exclude sterile neutrino scenarios. 
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